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Résumé 


Une analyse de la fonction ~Arg 6G + in) est présentée, cette analyse a permis de trouver une expression générale 


pour cette fonction qui n’utilise que des fonctions élémentaires comme la partie entiére et fractionnaire. Cette formule 


générale permet de mettre en lumiére une remarque générale qui avait été faite par Freeman Dyson a propos des 
E “7 : i A 1 in. ; . 
valeurs de la fonction ¢ reliées aux quasi-cristaux. On retrouve les mémes valeurs pour — Arg P'(- + > qui sont trés 
Tw 


voisines de la fonction ¢. De plus, les termes apparaissant dans ces formules générales obéissent a une régle simple 
avec n, a savoir, le nombre de 0 en binaire de n a partir de la droite. Ces approximations permettent également 


d’obtenir une formule plus exacte qui est asymptotique. Tous les résultats ont été obtenus de facon empirique. 


Abstract 


: : 1 ie ne. 8 : . . . ; 
An analysis of the function ~Arg ¢G + in) is presented. This analysis permits to find a general expression for that 
function using elementary functions of floor and fractional part. These formulas bring light to a remark from Freeman 
Dyson which relates the values of the ¢function to quasi-crystals. We find these same values for another function 
which is very similar, namely ~Arg rE+ >): These 2 sets of formula have a definite pattern, the n’th term is related 


to values like z,In(z) ,In(2), ..., log(p), where p is a prime number. The coefficients are closed related to a certain 
sequence of numbers which counts the number of 0’s from the right in the binary representation of n. These 
approximations are regular enough to deduce an asymptotic and precise formula. All results presented here are 


empirical. 
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Introduction 


Riemann établissait une formule qui compte le nombre de zéros de la fonction ¢ qui 
lorsque T tend vers I’infini est 


T T 7 


L’argument étant que si T est sur la droite % ou non, a été repris par von Mangoldt, Edwards et 
Hardy-Littlewood. Nous n’essayerons pas du tout d’en établir la preuve mais simplement de 
savoir combien de zéros y a-t-il dans chaque intervalle unité. En faisant la différentielle, en 
principe on obtient le nombre de zéros, on pose donc T > T + 1 et on fait la différence pour 
s’apercevoir que le terme d’erreur est trés imprécis. Tellement qu’en fait pour savoir combien y 
a-t-il de zéros dans chaque intervalle, la formule est inutile. 


Un certain Backlund en 1912 eut l’idée d’établir plus précisément que 
1 
NT) =~ HT) +14 S(T), (2) 


ou 0(T) est la fonction de Riemann-Siegel et S(T) est < Arg 65+ iT). Restait a trouver comment 


varie cette fonction d’angle avec la fonction ¢. Mais il se trouve que cette fonction S(T) varie 
énormément et a donné lieu a plusieurs idées de la part de Freeman Dyson, Alain Connes, Hugh 
Montgomery et M.V. Berry. (voir la bibliographie). Et plus récemment a André LeClair et 
Guilherme Franga de Cornell. 


Définissons F(n), comme étant le nombre de zéros de la fonction ¢ dans chaque intervalle unité. 
Ce qui nous donne la suite pour n = 1, 


F(n) = 0,0,0, 0,0, 0,0, 0, 0, 0,0, 0,0, 1, 0,0, 0,0, 0,0, 1,0,0,0,1,0,..., (3) 


la valeur est 2 a n=111, 150, 169, 223, et est 3 a n=5826, 5978, 6494. Il y a donc 2 zéros entre 
111 et 112. On compte ici le nombre de valeurs dont la partie entiére est 111. Comme on le sait 


: : 27M 
le nombre de zéros est environ : 
log(n) 
log(n) 


F(n) devrait étre de oe toujours lorsque n> oo. La valeur étant trés lentement croissante, il 


est aisé d’utiliser un dessin représentant les valeurs en tons de gris, 0 = blanc, 1= gris clair, etc. 
Etant donné qu’avec 1 image on peut représenter des milliards de valeurs. 


Ce qui équivaut a dire que pour chaque entier la valeur de 
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Ce dessin saute aux yeux, on y voit un motif immédiatement. Si on regarde attentivement le 
début sur la premiére ligne de pixels ont voit les zéros a 14, 21, 25, 30. La grille a été calibrée a 
4000 pixels de largeur, le choix est arbitraire et le motif s’y trouve a toutes les valeurs. Si on fait 
le dessin aux valeurs de n>> 1 on trouve (colorisée). 


Le motif de vague est clair, c’est un Moiré. Les zéros de la fonction ¢ sont des vagues. Ce n’est pas 
si surprenant, Granville l’a mentionné dans un article trés intéressant et aussi Baillie. Les zéros 
peuvent étre vus comme des ondes, en les sommant c’est une superposition d’ondes. 


En variant le pas de l’image entre 4000 et 9064 ou 9065, on peut percevoir le premier 


battement, la valeur 9064 ou 9065 est en fait le premier terme d’une superposition de quelque 
27 


log(2) 
termes ? M. V. Berry a suggéré log(p) ou des puissances de log(p), p étant premier. I] suffit de 
rajouter ces termes mais comment ? 


chose dont on ignore la nature, on a 9064 < 1000 < 9065. Mais quels sont les autres 


Le nombre de zéros d’une fonction 


Reprenons le calcul du nombre de zéros en commengant avec les fonctions les plus 
simples qui soit sont périodiques ou qui ont un grand nombre de zéros. La plus simple est sin(x), 
de facon évidente les zéros de sin(x) sont aux multiples de z, on prend ici x = 0. Si on construit 
la fonction comme F(n) alors on aura : 0,0,1,0,0,1,0,0,1, ..., en définissant 


h(n, a) = [(n + 1)a] — [na]. (4) 


: ‘ - ' 1 ‘anaes A . 
Ici [x] est la partie entiére de x. On a avec sin(x) que a = = Ce cas est trivial bien sir mais nous 


servira plus loin. Prenons un autre exemple avec la fonction de Bessel Jy) (x). Les zéros se situent 
a intervalles d’a peu prés 7, le premier est 2.40483, 5.52008, 8.65373, 11.7915, 14.9309, ... et 
pour ce qui nous intéresse, la suite 2,5,8,11,14, ... Sion applique donc le méme principe qu’avec 
sin(x) on trouve pas trop difficilement que le nombre de zéros B(n) de la fonction de Bessel J) (x) 
est donné par 
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B(n)=h (n, aa) (5) 


Ce qui est conforme a ce qu’on connait de cette fonction dont les zéros sont espacés d’un 
multiple de 7, petit a petit. Notons que la nature des zéros de cette fonction est inconnue. Ils ne 
semblent pas s’exprimer avec des constantes connues. 


Prenons un exemple plus difficile comme la fonction de Airy, nous prendrons la valeur positive 
pour simplifier l’argument, appelons A(n) la fonction qui compte le nombre de zéros par 
intervalle unité. Alors, 


2n'/2 
A(n) =h ( 37 (6) 


Celle-ci a été obtenue en utilisant la formule asymptotique connue (voir DLMF) et en utilisant la 
méthode dite bootstrap. Méme chose ici, on ne connait pas la nature des zéros de la fonction 
d’Airy mais on peut les compter par intervalle. 


La question est maintenant, peut-on trouver une expression semblable avec ¢ ?, si on utilise la 
formule asymptotique ¢a ne fonctionne pas, c’est trop imprécis. L’auteur a essayé une variété de 
méthodes bootstrap avec le log des premiers sans succés. Méme si certaines pistes semblaient 
prometteuses. Si y, est la partie imaginaire du n’iéme zéro non-trivial de la fonction ¢ alors 


2 Z 3 z 
¥1 (e™ —1)’ y2 (e” — 1) (7) 
Mais la s’arréte le motif, ce n’est qu’une heureuse coincidence. 


Densité ponctuelle 


On appelle densité ponctuelle, le nombre de valeurs d’une fonction dans un intervalle 
donné, ici l’intervalle est l’unité. Pour obtenir la densité ponctuelle ou une fonction qui 
représente la densité ponctuelle d’une autre il est nécessaire d’inverser fonctionnellement. 


n 
log(n) 
On peut bien sir étre plus précis comme le mentionne Dusart dans 


Exemple, le nombre de nombres premiers est et inversement, la densité ponctuelle des 


nombres premiers est : 
log(n) 


. : Sse 2 
sa thése. Pour la fonction ¢@, la densité est m 


log(ny’ 
Quand le terme principal est simple, la formule est directe mais si la formule est plus élaborée, 
c’est plus complexe. En reprenant la formule citée en [AL-GF] formule (13) on a le n’iéme zéro 
y(n) de la fonction ¢ est la solution de |’équation transcendante. 


; ae l 
ce qui donne une densité ponctuelle de oe) 


ym), ym) 1 ier sagt ah 
O77 bare +7 4rg CG + iv) es (8) 


Si on résous, on obtient directement avec Maple ou Mathematica 
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(n—-% 


=~ 21 ——_—_—_. 
y(n) - (@ = 2) (9) 


e 


La fonction W est la dite ‘de Lambert W’, quia des noms divers comme ‘ProductLog’. La formule 
est étonnante et trés précise quand n est grand, LeClair et Franca donnent des exemples ou n est 
de l’ordre de 101° et 107°° avec une précision qui dépasse de loin tout ce quia été fait avant. Ils 
notent également qu’ils ne tiennent pas compte des points de Gram ni non plus de la formule de 
Riemann-Siegel. Mais pour ce faire ils utilisent des méthodes numériques inversées, cette 
formule est le terme principal. La partie qui reste serait l’inverse fonctionnel du 
terme = Arg ¢(¢+in). Lerreur avec la formule (9) seule est comprise entre -1 et 1. Les 
méthodes utilisées précédemment ne permettent pas de dégager une formule directe pour 
corriger le tir. 


1 a 
Expressions pour — Arg ¢ (4+ in) 


Si on revient a la formule citée avec l’angle de la fonction “Arg ¢ G + in),en faisant le graphe on 


obtient, entre 1 et 115, 


-1 + 


Les valeurs sont assez imprévisibles. De plus, l’argument de la somme de deux nombres 
complexes méme simples se complique assez vite. Par exemple, en posant n=4000 et en faisant 


les somme partielles de “Arg f(¢+ i4000), la somme a l’infini est -0.382343520341 mais le 
graphe des sommes partielles nous donne: 
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HH! hi | LIP eH APA HAA Mt og 


4 Hit POPUL AT ! 


2 exemples de sommes partielles de <Arg c(2 + in) 


Donc, méme si la somme ultimement est -0.382343520341, pour s’y rendre les valeurs bougent 
bien trop pour prédire quoi que ce soit. Cette courbe n’est pas caractéristique du comportement 
général, d’autres valeurs de n donnent des courbes complétement différentes. 


On revient a la fonction qui donne approximativement le nombre de zéros, 


T TT 7 
NGS ae 10 
(T) ~ sr loss tg um) 


Si on rajoute le terme ~ Arg ¢(; + in), ou n remplace le T habituel, il apparait évident que la 


somme des deux termes tend vers une valeur qui est toujours prés de % en prenant la partie 
fractionnaire et ce de facon uniforme. Ce résultat est mentionné par LeClair et Franga. Mais, en 
regardant de prés, la valeur ajoutée (l’angle avec ¢), si on pose que c’est exactement ¥% alors ¢a 
implique que pour n=1 par exemple on a en (11) la valeur (qu’on appellera g(n)), une fois 
simplifiée nous donne, 


i 
g(1) = — 5 (4In(2) + 4 In(m) + 4 — 110) = 0.92333784.. (11) 
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La valeur de * Arg 6+ i) donne -0.43737200, la somme est 0.48596584 qui est proche de %, 


donc la valeur approchée de * Arg CG + i) serait = (4 In(2) + 4Iln(m) + 4 — 711). 


1 
g(2) = ieee In(t) + 8 — 11m) = 0.69231130... (12) 
En appliquant le méme principe on trouve donc la valeur approchée de “Arg oC +2i) = 
—0.19231128 qui comparée a la valeur exacte -0.19597758. Pour n = 10000 on trouve la valeur 
approchée de * Arg 6G + 10000i)= 


1 
— g— (120000 log(2) — 160000 log(S) + 40000 In(r) + 40000 + 81129 m) (13) 


On remarque que les valeurs sont des combinaisons linéaires de 1, log(t), log(2) et plus 
généralement log(p), avec p premier. On remarque également qu’a chaque valeur on se 
rapproche de la valeur supposée de ¥%. L’erreur a n=10000 est 0.0000066314559660306, cette 
valeur est - en fait, multiplié par un facteur de 24. Plus précisément, si e(10000) est l’erreur a 


n=10000, alors 


e(10000) = 0.0000066314559660306549805347837027 ... et 


240000 = 0.0000066314559621623056570368234738... (14) 


L’erreur est alors de l’ordre de 10~*° il reste suffisamment d’informations pour deviner le 
prochain terme qui est 7/120. Si on revient a la valeur trouvée en n=1, 2,3, on voit que cette 
correction est valide a partir de n=2, a n=10000000, l’erreur avec les 2 termes correctifs est de 
l’ordre de 10729. Le 3™e terme du développement asymptotique est 93/7!, le 4¢me est 127/896. 


Finalement, en reprenant chaque coefficient trouvé on peut écrire, 


8 


 ) | “(z 7 31 127 


1 1 n 
—Argé(2+i = {1 le (ap to tae t «). 
Tt res (5 in) nt Obes C\aen * 5760n3 * B0640n5 * 430080n7 )- (45) 


Et l’erreur est trés petite par rapport a n, pour n=10* ona une erreur de l’ordre de 10781. 


On reconnait ici le développement asymptotique de la formule de Riemann-Siegel, toujours avec 
la fonction {} qui représente la partie fractionnaire. On peut consulter les 2 suites A036282 et 
A114721 de la référence OEIS a ce sujet. 


: : “ges rar : 
On peut noter simplement également que les différences secondes de la suite —Arg 46 + in) est 


soit -1, 0 1 si on arrondit et c’est en analysant l’erreur par rapport a ces valeurs exactes qu’on 
peut utiliser la méthode de bootstrap pour trouver les coefficients de la formule 15. 
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Voici un tableau des valeurs de * Arg 4E + in) pour n=1-19. 


“Arg ¢ G A in) vraie valeur approximation valeur app. 
1 -0.437372012317 = (4In(2) + 4In(n) + 4-77) -0.423337836994 
2 -0.195977582921 = (8 In(m) + 8 — 7m) -0.192311274139 
3 -0.046800452442 = (-12 In(3) + 12 In(2) + 12 In(m) + 12 — 7 7) -0.0445622398292 
4 0.0474416622590 = (-16 In(2) + 16In(m) + 16 — 7m) 0.0491062514140 
5 0.101231367921 = (-20 In(5) + 20 In(2) + 20 In(m) + 20 — 7m) 0.102560818219 
6 0.122861669195 =-(-24 In(3) + 24 n(n) + 24 — 7m) 0.123968719880 
7 0.117778121706 ~(-4 In(7) + 4In(2)+4In() + 4-7) 0118726607797 
8 0.0898404109029 = (—64 In(2) + 32In(m) + 32 — 7m) 0.090670102212 
9 0.0419297327910 —(-72 In(3) + 36 In(2) — 36 n(@) + 36-77) 0.042667093960 
10 -0.0237198979997 —(-40 In(S5) + 40 In(t) + 40 — 771) -0.023056364340 
11 -.105325100472 = (—44 In(11) + 44 In(2) — 44 n(n) + 44-7 1) -0.104721949447 
12 -.201429006842 —(-48 In(3) — 48 In(2) + 48 In(n) + 48 — 70) -0.200876161160 
13 -.310818966587 aa (82 In(13) + 52 In(2) + 52 In(m) + 42 —7 1) -0.310308678190 
14 -.432469802098 a (-8 In(7) + 8 In(t) + 8-1) -0.431995985476 
15 434496598552 — (—60 In(3) — 60 In(S) + 60 In(2) + 60in(m) + 60 | 9 434938810386 

+10) 
16 0.290840842657 — (-192 In(2) + 64In(m) + 64+ 7) 0.291255403206 
17 0.137228161449 —(-68 In(17) + 68 In(2) + 68 In(m) + 68 + 1) 0.137618325910 
18 -0.0257546940666 = (-144 In(3) + 72 n(n) +72 + W) -0.025386213458 
19 -0.197586328497 =(-76 In(19) + 76 In(2) + 76 In(n) + 76+) -0.197237248073 


Expressions pour “Arg ¢ (3 - in), partie 2 


Les coefficients apparaissant dans la 3&me colonne sont particuliers. On a toujours le terme en 


87, vient ensuite le terme en mT, log(m),log(2) et log(p), les premiers p, sont ceux de la 


factorisation de n. 
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coefficient formule remarque 
log(1) 4n 
1 4nlog(2) + 4nlog(m) — 4nlog(n) + 4n — 71 
Tt 8 | = | - 77 
8 T 
4n sin =1mod 4 
—a(8 
log(2) a(n) = Osin =2 mod 4 a(8n) _ 4085058 
4n sin =3mod 4 32n 
_ (Osin=10u2 mod 3 -a(3n) _ 
log(3) a(n) = { pe eeey S) =A051064 


Le coefficient de log(2) est surprenant, si on fait le graphe des valeurs on obtient, 
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Les pics négatifs sont aux puissances de 2, les valeurs viennent en groupe de 3, suivis d’un zéro, 
la suite est 4, 0, 12, -16, 20, 0, 28, -64, 36, 0, 44, -48, 52, 0, 60, -192, 68, 0, 76, -80, 84, 0, 92, ... a(n) 
est parsemée de trous et a un maximum aux puissances de 2. En utilisant GFUN on trouve le lien 
assez facilement avec la suite A001511 ou A007814 du livre ou site de l’OEIS (voir 
bibliographie). La suite A001511 compte le nombre de zéros a partir de la droite en binaire de 


—-a 


<& ona la suite 2, 3, 2, 4, 2, 3, 2, 5, 2, 3, 2, 4, 2, 3, 2, 6, qui est A085058 du 


catalogue OEIS, c’est une variante de la suite générique A001511. 


n. Si on prend 


Le coefficient de log(3) suit la méme régle en décalé, la suite donne: 0, 0, -12, 0, 0, -24, 0, 0, -72, 
0, 0, -48, 0, 0, -60, 0, 0, -144, 0, 0, -84, 0, 0, -96, 0, 0, -324, 0, 0, -120, 0, 0, -132,... qui 
graphiquement donne, 
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C’est le méme motif fractal des suites A001511 ou A007814. Les autres termes en log(5), log(7) 


et plus généralement log(p) suivent le méme motif exactement mais de plus en plus décalés. Si 


on prend un terme sur 3 et on pose ——— — ar 


1, 1, 2, 1,1, 3,1, 1, 2, 1, 1, 2, 1, 1, 3,1, Fi ai est la suite A051064 qui est de la méme famille que la 
suite précédente. Les termes des autres rangs comme a(3n-1) et a(3n-2) sont triviaux ou 0. Le 
cas des coefficients de log(5), on retrouve la suite A091799 du catalogue OEIS, c’est une 


inspiré par les résultats précédents on trouve: 1, 1, 2, 


conjecture non vérifiée. 


On peut construire le méme tableau pour la fonction — ~ Arg rG +: it) qui suit les mémes régles. En 


fait, les 2 graphiques superposés de la fonction d’ cup de ¢ et del sont trés proches, l’une 
complétant l’autre, dans [AL-GF], on donne la fonction équivalente a g(n) serait 


ip) _ nlog Vm) 
a ae (16) 


Ces 2 fonctions sont les ondes porteuses de ce phénomeéne + la partie avec l’argument de ¢ et de T 
qui suivent les régles énoncées plus haut. 


i 1 1 1 in 
Approximations de — Arg r+: 


i ; 
Arg r (; + *) Valeur Approximation Valeur 


—1 
1 -0.380438567847 an (4ln(2) + m+ 4) -0.394472743168 
TT 


Etant donné que la fonction g(n) + ~Arg 6G + in) compte le nombre de zéros, si on ne considére 


que la partie fractionnaire en Arg c(t + in) c’est donc dire qu’a un multiple de m prés ¢a compte 


10 
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: : 1 : ’ : 
aussi le nombre de zéros. Comme on le sait, — Arg c(2 + in) varie entre -1 et 1. Donc, si effectue 


les différences, A de cette fonction on obtient le décompte aussi, méme si la valeur exacte 
comme tel est pres d’un entier +1/2, il en est de méme pour A., a partir de n=1, g(n) + 


“Arg c(2 + in) - 


0.485966, 0.496333, 0.497764, 0.498338, 0.498669, 0.498892, 0.499051, 0.499168, 0.499263, 
0.499335, 0.499395, 0.499447, 0.499489, 0.499526, 1.49956, 1.49958, 1.49961, 1.49964, 
1.49965, 1.49967, 1.49968, 2.49970, 2.49972, 2.49972, 2.49974, 3.49974, ... 


Qui arrondis et en faisant la difference donne bien le compte des zéros de la fonction F(n) 
originale 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, ... qui comme le mentionne 
[AL-GF] manque une valeur entre 1007 et 1008 parce que + ~ Arg cE + in) fait 1 tour de plus. La 


fonction donne bien |’angle exact mais au multiple de mt prés. Il en est de méme avec la valeur 
approchée. 


Dans les 2 cas, ¢a décompte bien le nombre de zéros, si la valeur différe, il y a un tour ou 
plusieurs qui sont effectués par l’une ou I’autre expression mais le fond reste le méme. Ce qui a 


été montré ici est que le grain ou la variation mystérieuse de la fonction —Arg ¢(5 + in) et 

1 i : celccah 

Arg ré + in) sont en fait des approximations de constantes comme T, log(t), log(2) et log(p) 

avec l’aide de la fonction [ ] et { }. Donc, pour les premiéres valeurs de n et pour beaucoup 
; ; j : pee! n ae . 

d’autres, on a bien le décompte des zéros mais, a mesure que n croft, la fonction Arg 4E + in) 

fait de plus en plus de tours a chaque entier. 


Ce qui rejoint assez exactement les propos de Freeman Dyson a propos de quasi-cristaux. I] 


suggérait qu’a la base, les zéros ou les valeurs de la fonction ¢ sont des variations générées a 
l’aide de constantes jumelées avec ces 2 fonctions [ ] et { }, que ce que l’on voit est un reflet. 


11 


é 1 . 1 in 
On the values of the functions — Arg ¢(t+in) and —Arg r+ 


Bibliographie et sur internet 


[AMO] A. M. Odlyzko, On the Distribution of Spacings Between Zeros of the Zeta Function, 
Mathematics of Computation 177 (1987) 273. 


[HME] H. M. Edwards, Riemann's Zeta Function, Dover Publications Inc. 1974. 
[ECT] E. C. Titchmarsh, The Theory of the Riemann Zeta-Function, Oxford University Press, 1988. 


[MVB] M. V. Berry, Riemann's Zeta Function: A Model for Quantum Chaos?, Quantum Chaos and 
Statistical Nuclear Physics, Eds. T. H. Seligman and H. Nishioka, Lecture Notes in Physics, 263 
Springer Verlag, New York, 1986. 


[AL-GF] André LeClair and Guilherme Franga, Statistical and other properties of Riemann zeros 
based on an explicit equation for the n-th zero on the critical line. Arxiv [1307.8395] et la version 
2 également Arxiv [1307.8395v2]. 


[XG] Xavier Gourdon, The 10?3 first zeros of the Riemann Zeta function, and zeros computation at 
very large height, (2004). 


[AG] An electrostatic depiction of the validity of the Riemann Hypothesis and a formula for the N-th 
zero at large N, article original sur Arxiv : arXiv:1305.2613. 


[AG] Andrew Granville, Prime numbers patterns, 
http: //www.maa.org/sites /default/files /images /upload library/22/Ford/Granville.pdf 


[AG2] Andrew Granville, Prime numbers race, 
http: //www.dms.umontreal.ca/~andrew/PDF/PrimeRace.pdf 


[DLMF] La fonction d’Airy http://dlmfnist.gov/9.7 


[FD] Freeman Dyson, http://empslocal.ex.ac.uk/people/staff/mrwatkin/zeta/dyson- 
quasicrystal.htm 


[FD2] Freeman Dyson, Birds and Frogs, 2009 : 
http://www.ams.org/notices/200902 /rtx090200212p.pdf 


[OEIS] Neil James Alexander Sloane The On-Line Encyclopedia of Integer Sequences 


[EIS] NJA Sloane et Simon Plouffe, The encyclopedia of integer sequences, Academic Press, 1995, 
pp 600. 


[BP] Francois Bergeron et Simon Plouffe, GFUN, F. Bergeron and S. Plouffe, Computing the 
Generating Function of a Series Given its First Few Terms, Experimental Mathematics, Volume 1, 


(1992), 307-312. (ps.gz, pdf ). 


[SP] Approximations of generating functions and a few conjectures, arxiv 
http: //arxiv.org/abs/0911.4975 


[WB] La fonction de Bessel http: //fr.wikipedia.org/wiki/Fonction de Bessel 


12 


é 1 . 1 in 
On the values of the functions — Arg ¢(t+in) and —Arg r+ 


[WZ] La fonction Zeta de Riemann : 
http: //fr.wikipedia.org/wiki/Fonction z2%C3%AAta de Riemann 


[WA] La function d’Airy http://fr.wikipedia.org/wiki/Fonction d%27Airy 


[LMFDB] Les 103 800 788 359 premiers zéros calculés par David Platt 
http: //www.|mfdb.org/zeros/zeta 


[PDU] Thése de Pierre Dusart, http://www.unilim.fr/laco/theses/1998/T1998 01.html et 
http://www.unilim.fr/laco/theses/1998/T1998 01.pdf 


[RB] demonstration project on the zeta function by Robert Baillie : How the Zeros of the Zeta 
Function Predict the Distribution of Primes. 


[RB] Robert Baillie, 
http://demonstrations.wolfram.com/HowTheZerosOfTheZetaFunctionPredictTheDistributionO 


fPrimes/ 


[RB] Robert Baillie, 
http://demonstrations.wolfram.com/author.html]?author=Robert%20Baillie 


[AS] J.-P. Allouche and J. Shallit, The ring of k-regular sequences, Theoret. Comput. Sci. 98 (1992), 
163-197. 


[AS2] J.-P. Allouche and J. Shallit, The Ring of k-regular Sequences, IJ, Theoret. Comput. Sci. 307 
(2003), 3-29. 


[EW] The Riemann-Siegel Theta functions http://mathworld.wolfram.com/Riemann- 
SiegelFunctions.html 


[LS] J. C. Lagarias and N. J. A. Sloane, Approximate squaring , Experimental Math., 13 (2004), 113- 
128. 


13 


